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Analysis of the low-energy electron-recoil spectrum of the CDMS experiment 
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We report on the analysis of the low-energy electron-recoil spectrum from the CDMS II experiment 
using data with an exposure of 443.2 kg-days. The analysis provides details on the observed counting 
rate and possible background sources in the energy range of 2 - 8.5 keV. We find no significant excess 
in the counting rate above background, and compare this observation to the recent DAMA results. 
In the framework of a conversion of a dark matter particle into electromagnetic energy, our 90% 
confidence level upper limit of 0.246 events/kg/day at 3.15 keV is lower than the total rate above 
background observed by DAMA by 8.9cr. In absence of any specific particle physics model to provide 
the scaling in cross section between Nal and Ge, we assume a Z 2 scaling. With this assumption 
the observed rate in DAMA differs from the upper limit in CDMS by 6.8<r. Under the conservative 
assumption that the modulation amplitude is 6% of the total rate we obtain upper limits on the 
modulation amplitude a factor of ~2 less than observed by DAMA, constraining some possible 
interpretations of this modulation. 
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PACS numbers: 29.40.-n, 95.35. +d, 95.30.Cq, 85.25.Oj, 29.40.Wk 



Astrophysical observations strongly suggest that non- 
luminous, non-baryonic matter constitutes most of the 
matter in the Universe. This dark matter should be lo- 
cally distributed in dark halos of galaxies such as the 
Milky Way, enabling the direct detection of the dark mat- 
ter particles via their interaction in terrestrial detectors. 
The movement of the Earth around the Sun would pro- 
vide an annual modulation of the counting rate, caused 
by the change in the relative velocity of the dark mat- 
ter particles and the earthbound target. The DAMA 
collaboration claims the observation of such a modula- 
tion, in two different Nal(Tl) scintillation detector ar- 
rays, the original DAMA/Nal setup [l| and the upgraded 
DAM A/LIBRA experiment [2j. The observed signal is in 



the 2-6 keV electron-equivalent energy range with a peri- 
odicity of 0.998±0.003 years and a phase of 144±8 days. 
The DAMA collaboration claims that no known system- 
atic detector effect could explain the modulation signal. 
The modulation phase is consistent with the expected sig- 
nature of galactic dark matter particles interacting in a 
terrestrial detector. However, the original interpretation 
of the DAMA result as a signal from Weakly Interact- 
ing Massive Particles (WIMPs) that would interact via 
nuclear recoils [l| is inconsistent with other experimental 
results 0,11 HE 0| Hi- Note, that the DAMA detectors 
do not discriminate between electron recoils and nuclear 
recoils. 

A signal from an electromagnetic dark matter interac- 
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tion should be detectable in the Cryogenic Dark Matter 
Search (CDMS) experiment, but would be rejected in our 
standard search for nuclear recoils [|| . The possibility of 
an electron-recoil signal from axion-like dark matter par- 
ticles has recently been investigated || @, In this 
paper, we present a general analysis of our low-energy 
electron-recoil spectrum, provide details on the observed 
counting rate in this energy range, and comment on the 
implications of these results on possible interpretations of 
the energy spectrum and the modulation signal observed 
by DAMA. 

The CDMS collaboration operates a total of 19 Ge and 
11 Si crystal detectors, each having a mass of ~250g and 
~100g, respectively, at a temperature of ~40mK in the 
Soudan Underground Laboratory 11,13]. The ionization 



and phonon energy of every event is read out simultane- 
ously. The recoil energy is reconstructed from them. The 
ratio of ionization to recoil energy, the ionization yield, 
discriminates nuclear- from electron-recoil events. 

In this analysis we consider data with a total exposure 
of 443.2 kg-days before cuts, which has been acquired 
in two run periods between October 2006 and July 2007 
(designated as R123 and R124) and is the same dataset 
used for an axion search analysis [lfj. Three of the 19 Ge 
detectors were excluded because of readout failures and 
another one due to reduced trigger performance at low 
energies. From the remaining 15 Ge detectors one suf- 
fered from reduced trigger performance in R123 and two 
from incomplete neutralization in R124 which have also 
been left out of the analysis. The silicon detectors were 
not considered. We required that an event had to pass 
several cuts. The events needed to have an ionization 
energy at least 3a above the mean noise and be recorded 
in only one detector. All 30 detectors were used to se- 
lect these single-scatter events. Moreover, we demanded 
that there was no signal in the scintillator veto shield 
surrounding the detectors. The length of the veto coin- 
cidence window was set to 50 /is. In order to explore the 
low-energy electron-recoil spectrum we selected events in- 
side the 2(7 electron- recoil band in ionization yield 
The fiducial volume was measured using nuclear-recoil 
events from calibrations with a 252 Cf source because of 
the uniform distribution of neutrons throughout the de- 
tector. We excluded all datasets taken within 3 days af- 
ter a neutron calibration to avoid high gamma rates due 
to neutron activation of the detectors' supporting struc- 
ture made predominantly of copper. The remnant rate of 
64 Cu contributes less than 2% to the mean counting rate 
at low energies and decreases with a half life of 12.7 h. 

The summed background spectrum of all considered 
detectors, taking into account the detection efficiency 
[Toj ] . is shown in Fig. [I] For reference, the corresponding 
counting rates are also given in TablelH In this analysis 
we consider the electron-equivalent energy range between 
2 and 8.5 keV based on the ionization signal, in which 
the background rate is ~1.5 events/kg/day/keV. Fig.[T] 
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FIG. 1: Fit (red line) to the efficiency corrected low-energy 
spectrum consisting of a background model (gray/dashed) 
and three Gaussian distribution functions describing the 
10.36 keV line from 71 Ge (black), the 8.98 keV line from 65 Zn 
(blue) and a line at the energy of 55 Mn (green, see text). The 
total counting rate of the latter two lines is given in the figure. 

also illustrates a simple fit to the observed electron-recoil 
spectrum. The fit incorporates known spectral lines at 
10.36 keV and 8.98 keV, both outside of our analysis win- 
dow. The former is caused by X-rays and Auger-electrons 
from the decay of 71 Ge, a product of neutron capture on 
70 Ge during neutron calibrations. The latter originates in 
the decay of remnant 65 Zn from cosmogenic activation of 
the detectors. We also fit for a spectral line correspond- 
ing to an excess of events observed near 6.5 keV, which 
is likely caused by the de-excitation of 55 Mn; this feature 
is discussed further below. Each peak is fit by a Gaus- 
sian distribution function with width fixed at CDMS's 
measured energy resolution [lij. The detector-averaged 
r.m.s. energy resolution a(E) below 10 keV is given by: 



a{E) = v/(0.293) 2 + (0.056) 2 £ [keV] 



(1) 



where E is the measured energy in keV. 

55 Mn can be produced by electron capture of remnant 
55 Fe from cosmogenic activation. The de-excitation of 
55 Mn results in a spectral line at 6.54 keV, matching ex- 
actly the energy of the corresponding peak in our spec- 
trum. While at the surface the detectors were exposed 
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TABLE I: Rate [events/kg/day/keV] in the 2 - 8.5 keV energy 
range. 
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to fast neutrons from cosmic-ray showers. Gamma-rays 
from isotopes produced in Ge by these fast cosmic-ray 
neutrons have been observed in the CoGeNT experiment, 
which uses a p-type contact germanium detector provid- 
ing an excellent energy resolution Q. The most domi- 
nant lines in their spectrum are from 65 Zn with an en- 
ergy of 8.98 keV and 68 ' 71 Ge with an energy of 10.36 keV, 
which arc both also visible in our spectrum. Calculations 
of the production rate of cosmogenic isotopes show that 
55 Fe is produced in Ge 0. The long halflife of 55 Fe of 
2.73 y allows a remaining activity of this isotope in the 
detectors. Since the activation stopped when the detec- 
tors were moved underground, the time evolution of this 
counting rate would enable us to determine if it is caused 
by 55 Fe isotopes. However, the uncertainties in the pro- 
duction rate and on the time the detectors spent at the 
surface are likely too great to give a reliable constraint on 
the total rate expected from the de-excitation of 55 Mn. 

We carried out a profile likelihood analysis in order to 
search for an excess of event rate above background [3] • 
The event rate per unit measured energy (E) and per 
detector (d) including background was written as: 

R(E, d) = B(E, d) + A(E, d) (2) 

The background B(E, d) is assumed to be of the form 

H(dY 



B(E, d) = e{E,d) 



C{d) + D(d)E + 



7]-£(E,d) 



A, 



0.54 



E 



(3) 

/ 27r<7 6 . 54(d) 

where C(d), D(d) and H(d) are free parameters de- 
termined by the fit routine and e(E, d) is the energy- 
dependent detection efficiency. The Gaussian represents 
a contribution from 55 Fe decays at an energy of 6.54 keV. 
A(E,d) represents a spectral line at a given energy Eq. 
Thus, we used a Gaussian distribution function multi- 
plied with the efficiency: 

A 



A(E,d) = e(E,d) 



. V2<j (d) , 



(4) 



/27n7 (cO 

Since we have no constraint on the 55 Fe contribution to 
the spectrum we do not subtract a possible background 
contribution. The reason for introducing the additional 
factor 77 in ([3]) is that, while scanning over the recoil en- 
ergy and approaching the 6.54 keV background peak, the 
fit function actually consists of a sum of two Gaussians at 
the same energy. Thus, it serves as a weight suppressing 
the importance of the 55 Fe rate in the background model 
B(E,d). We varied 77 in steps of 0.1 between and 1 
and took the most conservative of these limits for each 
energy. 

The fit was performed by a maximization of the un- 
binned log-likelihood function 
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FIG. 2: Fit (red) to the published DAMA/LIBRA low-energy 
spectrum @], consisting of a background model (grey/dashed) 
and a Gaussian distribution function (green). The parameters 
of the Gaussian are given in the figure. 

where the sum goes over events (i) and detectors (J), with 
respect to Ao and the background parameters. Rt de- 
notes the total sum of the event rate (R) over energy and 
all detectors. We find no statistically significant excess of 
the event rate above background. We set a Bayesian 90% 
confidence level (CL) upper limit on the total counting 
rate A by integrating the profile likelihood function in 
the physically allowed region (A > 0). 

The annual modulation signature observed by DAMA 
Q may be interpreted as the conversion of a dark matter 
particle into electromagnetic energy in the detector. In 
this case it should be possible to observe the correspond- 
ing signal in the electron-recoil spectrum of CDMS. The 
upper limits on an excess rate presented in this paper 
should thus help to identify or constrain possible mod- 
els which can explain the annual modulation signature 
observed by DAMA. The total counting rate above back- 
ground observed by DAMA/LIBRA in the claimed signal 
region has been obtained from a fit to their spectrum con- 
sisting of a Gaussian and a background model shown in 
Fig. [2]giving a total rate of 0.698 ± 0.051 events/kg/day. 
A direct comparison between the 90% CL upper limits 
from this analysis (black/solid) and the total rate above 
background observed by DAMA (black data point with 
2a error bars in the figure) is shown in Fig. [3] At the 
energy of the DAMA peak (3.15keV) the observed rate 
is 8.9<7 away from the upper limit on the rate in CDMS 
of 0.246 events/kg/day. Though the peak of Fig. O may 
contain a contribution from the decay of 40 K and the sub- 
sequent de-excitation of 40 Ar resulting in a spectral line 
at 3.2 keV, no information is supplied on the actual rate 
of such a background [l5j |. Thus, no subtraction is per- 
formed, which would reduce the difference between the 
upper limit from CDMS and the excess rate in DAMA. 

The event rates in the CDMS and DAMA detection 
media may differ depending on the coupling of the dark 
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FIG. 3: 90% CL upper limit on the total counting rate in 
Ge from this analysis (black/solid). The corresponding upper 
limit on the total counting rate in Nal under the assump- 
tion of a Z 2 scaling of the conversion cross section (see text) 
is also shown (blue/dashed). The black data point with 2a 
error bars gives the total counting rate of the 3.15 keV peak 
of DAMA/LIBRA derived from a fit to their spectrum (see 
Fig. [2}. The insert compares the upper limit on the mod- 
ulation amplitude assumed to be 6% of the unsealed up- 
per limit (black/solid) and the Z 2 scaled upper limit in Nal 
(blue/dashed) with the 2a regions of the annual modulation 
amplitude observed by DAMA (Nal+LIBRA) in the 2 - 4 keV 
(red/filled) and 2-6 keV (green/hatched) energy range. 

matter particle. Thus, the upper limits in Ge have to be 
scaled to the expected rate in Nal in order to perform 
a comparison in a particular model. For an electromag- 
netic conversion of a dark matter particle, the particle 
velocity is essentially irrelevant (in contrast to the calcu- 
lation for nuclear recoils, where the energy threshold pro- 
vides a minimum velocity for the phase space integral). 
Thus, the annual modulation signature is only caused by 
a change in the particle flux over the course of the year. 
The total counting rates per unit mass of such a conver- 
sion in the case of a Ge and a Nal target are related by 
the following condition: 

RNal _ Ace 07 + (Tjv Q 

Roe Aj + An 

where A4 is the atomic mass of the nuclei, and Ui is the 
total cross section per atom of the interaction. The de- 
tection efficiencies in both materials should be very close 
to 100% at these low energies; thus, effects of a material 
and detector geometry dependent detection efficiency are 
neglected in the following. 

The total cross section will depend on the coupling of 
the dark matter particle to the detection media. For an 
electromagnetic conversion a Z 2 (where Z is the atomic 
number) scaling of the cross section is natural and is thus 
considered in the comparison of the rate limits in Ge from 
this analysis with the rate observed by DAMA. Another 



scaling can be trivially considered. This is a more general 
model than the one considered in our axion search paper 
ficj |. The scaled rate limits in Nal at a 90% CL are 
given in Fig. [3] (blue/dashed line). The total counting 
rate observed by DAMA/LIBRA is 6.8c greater than the 
upper limit at 3.15keV. 

Under standard halo assumptions a conservative upper 
limit on the modulation amplitude is ±6% if the modu- 
lation is caused by a change in the particle flux only [l6| . 
Note, that if the conversion cross section is inversely pro- 
portional to the dark matter particle velocity (as inelas- 
tic cross sections tend to be 17] ) the annual modulation 
amplitude is highly suppressed. The insert in Fig. [3] com- 
pares the unsealed upper limit (black/solid) and the Z 2 
scaled upper limit in Nal (blue/dashed) on the modula- 
tion amplitude with the 2a regions of the annual mod- 
ulation amplitude observed by DAMA (Nal+LIBRA) in 
the 2-4 keV (red/filled) and 2-6 keV (green/hatched) 
energy range [2j . The upper limits on the modulation am- 
plitudes are a factor of ^2 less than observed by DAMA. 

In this paper we reported on our analysis of the low- 
energy electron-recoil spectrum of the CDMS experi- 
ment, providing the observed rate in the 2 - 8.5 keV range 
and the identification of possible background sources. 
The analysis found no significant excess in the counting 
rate above background. Considering the conversion of a 
dark matter particle into electromagnetic energy the 90% 
CL upper limit on the total counting rate from CDMS 
at 3.15 keV is 8.9<r (6.8<r) below the excess rate observed 
by DAMA in a direct comparison (under the assumption 
of a Z 2 scaling of the cross section), neglecting a pos- 
sible background contribution from 40 K in the DAMA 
data. We note that the actual scaling between Ge and 
Nal has to be provided by a specific model, but stress 
that an analysis of the low-energy electron-recoil spec- 
trum of CDMS helps to identify or constrain possible 
models which can explain the annual modulation signa- 
ture observed by DAMA. In the conservative case of a 
6% modulation amplitude our recent data provides 90% 
CL upper limits on the modulation amplitude that are a 
factor of ~2 less than observed by DAMA. 
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